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Abstract: In the present note, a summary of selected aspects of time- 

dependent mean-field theory is first recalled. This approach is optimized to 
describe one-body degrees of freedom. A special focus is made on how this mi¬ 
croscopic theory can be reduced to a macroscopic dynamic for a selected set of 
collective variables. Important physical phenomena like adiabaticity/diabaticity, 
one-body dissipation or memory effect are discussed. Special aspects related 
to the use of a time-dependent density functional instead of a time-dependent 
Hartree-Fock theory based on a bare hamiltonian are underlined. The absence 
of proper description of complex internal correlations however strongly impacts 
the predictive power of mean-field. A brief overview of theories going beyond 
the independent particles/quasi-particles theory is given. Then, a special atten¬ 
tion is paid for finite fermionic systems at low internal excitation. In that case, 
quantum fluctuations in collective space that are poorly treated at the mean- 
field level, are important. Several approaches going beyond mean-field, that are 
anticipated to improve the description of quantum fluctuations, are discussed: 
the Balian-Veneroni variational principle, the Time-Dependent Random Phase 
Approximation and the recently proposed Stochastic Mean-Field theory. Re¬ 
lations between these theories are underlined as well as their advantages and 
shortcomings. 
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1 Introduction 

The time-dependent mean-field approach provides a unique tool to describe nuclear dy¬ 
namic from small to large amplitude collective motion. Microscopic time-dependent Hartree- 
Fock (TDHF) based on zero range interaction [T] were introduced only four years after the 
seminal work of Vautherin and Brink [2] showing that Skyrme like interaction can be used 
to construct an energy functional of the density (Energy Density Eunctional [EDE]). The 
EDE leads to an accurate description of selected static properties in nuclei [am. Static 
and dynamical mean-field theories, by considering directly single-particle degrees of free¬ 
dom (DOE) interacting through a common mean-field, was a major breakthrough in nuclear 
physics and opened the perspective of a unique framework able to describe static, dynami¬ 
cal and thermodynamical nuclear properties. In particular it avoids the preselection of few 
collective DOE, generally necessary in macroscopic theories. Already at the very early stage 
of the application of Time-Dependent EDE, it was realized that the independent particle 
picture used in the theory leads to severe limitations. In particular, mean-field theory that is 
anticipated to be optimal for one-body DOE provides a rather poor description of two-body 
and higher-order degrees of freedom. Many aspects of nuclei, like shape evolution, relative 
distance of two nuclei during a collision, ... are associated to one-body collective observ¬ 
ables. Two-body or higher order DOE are connected to fluctuations in collective space. The 
poor treatment of quantum fluctuations leads to a quasi-classical description of collective 


motion within nuclear mean-field theory. This has important consequences on the predictive 
power of this approach. TD-EDF will be particularly powerful to describe specific physical 
phenomena where quantum fluctuations do not play a decisive role. It is for instance able 
to nicely predict the mean collective energy of giant resonances [5] . It also provides a very 
reasonable description of the fusion process above the Coulomb barrier. In recent years, we 
have seen an increasing interest in applying TD-EDF due to the possibility (i) to perform 
symmetry unrestricted evolution, i.e. full three-dimensional calculations (see other chapter 
of the present book) [H [3 0 [H [TOl [TTl [121 IE] (ii) due to the possibility to use complete 
Skyrme functionals compatible with state of the art nuclear structure calculations (iii) more 
recently, the possibility to incorporate pairing correlations musiiiiiiTiiiaiiaiio]. with 
these improvements, many applications ranging from collective vibrations, fusion, deep in¬ 
elastic collisions like multi-nucleon transfer or quasi-fission or more recently fission have 
shown that microscopic mean-field theories can provide important physical information. 
Unfortunately, in all cases, some shortcomings uncovering the poor treatment of quantum 
aspects in collective space recurrently shows up: 

• Life-time of giant resonances: the mean energy of giant resonances which are 
obtained using TD-EDF subject to small perturbation are usually properly described. 
This is usually not the case of its life-time. Indeed, collective motion are expected to be 
damped out due to the coupling of this excitations with more complex internal degrees 
of freedom. This coupling directly reflects correlations inside the nucleus and is not 
properly treated in TD-EDF. As a consequence, the spreading width of the nuclear 
response is usually severely underestimated (see for instance 0 ). Giant resonances 
associated to shape oscillations are associated to one-body degrees of freedom and 
therefore TD-EDF should a priori be an optimal approach for this process. While, as 
we will see, this is indeed the case for short time, other DOFs might affect the one- 
body dynamic for long time. The long-time description of nuclear systems requires to 
incorporate the effect of complex DOFs that has been disregarded at the mean-field 
level, at least in an approximate way. 

• Underestimation of experimental fluctuations: when TD-EDF was applied to 
deep inelastic collision, it was immediately realized that it can rather well describe the 
mean properties of the reaction partners like the average mass of the quasi-projectile 
or its average total kinetic energy but completely fails to describe the width of ex¬ 
perimental distributions. From the extensive comparison between experiments and 
theory, two important conclusions can be drawn. Since the average properties of 
collective variables are well describe and since important dissipation occurs in deep 
inelastic collisions, TD-EDF properly incorporate dissipation in one-body space. More 
precisely, it accounts for the so-called one-body dissipation effects that is associated 
to both the nucleon exchange process and to the reaction of single-particle states to 
a dynamical shape deformation of the nucleus. Experimental fluctuations are directly 
linked to fluctuations in collective space. These fluctuations that are associated to 
two-body operators, are beyond the scope of the theory. It should be noted that quan¬ 
tum fluctuations are not zero in a mean-field theory and therefore this approach might 
have given a reasonable account for the fluctuations. Unfortunately, nuclei are finite 
many-body quantum systems were large quantum fluctuations takes place and were a 
simple mean-field theory strongly underestimates them. 

• Tunneling in collective space: a last example of the limitation of TD-EDF is the 
complete absence of quantum tunneling in collective space. As a consequence, two 
nuclei cannot fuse if their center of mass energy is below the fusion barrier energy. 
Similarly, the spontaneous fission process that requires the fission barrier transmission 
cannot be described. 

These few examples illustrate the necessity in some cases to improve the microscopic 
description of small and large amplitude dynamics beyond the independent particle picture. 
In the last decades, several methods have been proposed to correct mean-field evolution. 


Most often, due to the increase of complexity compared to the original TD-EDF, these the¬ 
ories have been only applied to rather schematic models or even not been applied at all. 
Nevertheless, TD-EDF has now reached a certain level of maturity and with the constantly 
increasing computational power, the standardization of time-dependent methods beyond 
mean-field appears as the next natural step in the field. The aim of this article is to pro¬ 
vide some guidance in the beyond mean-field approaches dedicated to many-body fermionic 
motion. Several reviews have been recently made on this subject either focusing on the 
inclusion of correlations [5], or on the addition of new DOFs through a variational ansatz 
m or the use of stochastic methods [21]. The goal of the present article is not to redo an 
exhaustive review but rather to present selected theories that might be of practical interest 
in the coming years focusing especially on collective aspects. First, some general aspects of 
mean-field are recalled. Specific connections between single-particle and collective DOFs are 
introduced. Finally, several theories that are expected to improve mean-field and to treat 
quantum fluctuations in collective space are underlined. 

2 Mean-field approach and collective evolution 

The connection between mean-field and selection of one-body DOFs is best illustrated 
starting from a variational principle 

5 = ( 1 ) 

Jto 

where H is the hamiltonian while the variation of the action is restricted to Slater determi¬ 
nants. Due to the Thouless theorem [22], local variation in the Slater determinant Hilbert 
space can be written as: 

( 2 ) 

Said differently, the one-body DOF associated to the set of operators {a\aj} are the gen¬ 
erators of the local transformation between Slater determinants. This specific property in 
combination with the use of the variational principle 0 insures that one-body DOFs fulfill 
the Ehrenfest evolution ESIIII]: 

At this stage, many important comments can be made: 

• Equation 0 is valid all along the dynamical path. This does not mean that the 
one-body evolution is exact because there is no reason that the exact time-dependent 
state remains a Slater determinant (SD) at all time. Starting from a SD state, due to 
correlations, the exact many-body state will be a complex superposition of many SD. 
However, if the state is initially prepared as a SD, the one-body density match the 
exact evolution for short time. It is in that sense that the dynamic of one-body DOFs 
is optimized. 

• Eq. (|^ is a priori only valid for single-particle operators. In particular, the Ehrenfest 
equation is not fulfilled for the two-body density matrix components (dtatd/cd/) along 
the mean-field path. So even for short time, two-body and/or higher order effects 
might be poorly treated. 

• Since the expectation value in Eq. 0 is taken over a Slater determinant, due to the 
applicability of the Wick theorem, the right hand side of 0 becomes also a functional 
solely of the one-body density components. We then end-up with a closed set of 
equation that can schematically be written as: 

( 4 ) 




Denoting by pji = (nta^) the one-body density matrix component. The above equation 
identifies with the standard TDHF evolution: 

ihdtp = [h(p),p] (5) 

where h{p) is the self-consistent mean-field hamiltonian. 

The mean-field equation applies to more general situations where the system is not 
necessarily described by a Slater determinant but by a many-body density written as: 


D = exp I — 


ijdlcLj 


Zq 


(6) 


where Zq insure that the density matrix is properly normalized. A new variational 
principle should be introduced starting directly from a density matrix instead of a trial 
wave-function. Such a variational principle has been proposed by Balian and Veneroni 
in the 80’s [24l [25j [26] . For a recent comprehensive discussion see m- The connection 
between the dynamical path and the selection of few DOFs is directly visible in this 
variational principle. Given a set of observables, denoted generically by A = {A^} and 
denoting by A(t) and D{t) the observable and density matrix written in the Heisenberg 
representation. The action to minimize becomes: 


5 = Tr 




( 7 ) 


The action is varied by imposing the two boundary conditions D{ti) = Di^ A(tf) = 
A/. We see in particular from this variational principle, that the evolution will depend 
on the specific set of observables that are selected. In particular, when the action is 
varied on the specific class of densities given by Eq. © and the selected sets of relevant 
observables are one-body observables, i.e. A = Z[-a\aj^ the action becomes a 
functional of Z(t) and Z'(t), i.e. S = tS(Z(t), Z'(t)). the variation with respect to 
7a'( t) gives: 

= T,([d\a^,H\f){t)), (8) 


that match also the exact Ehrenfest dynamics expect that the density is now restricted 
to Due to the specific algebra of densities given by Eq. (|^, Wick theorem can 
also be applied and the mean-field Eq. © is recovered. This more general derivation 
applies in particular if the initial state is a statistical ensemble of independent particles. 
As we will see below, another interest of the BV variational principle is that it might 
also be used in cases where the observable of interest are not necessarily one-body 
observables. 


Pairing correlation: In the last ten years, efforts have been made to incorporate 
pairing correlations in mean-field applications. While technically more involved, the 
Time-Dependent Hartree-Eock Bogoliubov (TDHEB) theory is formally a straightfor¬ 
ward generalization of TDHE [4] where the independent particles picture is replaced 
by independent quasi-particles. Breaking the U{1) symmetry associated to particle 
number conservation leads to non-zero value of the anomalous density components, 
i^ij = {ajai). TDHEB approach is optimal for a generalization of one-body ’’pseudo¬ 
observables” that can now be written as (using the notation of [3|): 

i = ^Zl^aldj + ^[zf^aja, + Z]°*ala]). (9) 

ij ij 


Breaking the particle number symmetry, is an elegant way to include part of the 
correlations, i.e. long-range pairing correlation, in the description. Still it is often 






classified as a mean-field theory where the one-body density p and mean-field h[p] are 
respectively replaced by the generalized density IZ and generalized mean-field 
m- In particular, including pairing does not solve any of the shortcoming listed in the 
introduction. The recent systematic applications of refs. EZIIIH] clearly demonstrate 
the absence of damping in giant resonance. In addition, although TDHFB considers 
an enlarged set of relevant observables, their dynamics is still very classical. In the 
following, mean-field (MF) term will be used for approaches like TDHF or TDHF, 
while beyond mean-field (BMF) is reserved to approaches beyond the independent 
particles or quasi-particles. 

• Nuclear mean-field and density functional theory: Last, I would like to mention 
an important aspect of the nuclear mean-field that has sometime been underestimated. 
In the present section, mean-field is introduced starting from a well-defined hamilto- 
nian H. This situation will be referred here as the ’’Hamiltonian case”. Most often, the 
hamiltonian case can only serve as a guidance in the nuclear mean-field context. In¬ 
deed, what we call a mean-field approach is actually a density functional theory (DFT) 
where the energy is written as a functional of the one-body density matrix (normal 
or generalized) and where the parameters of the functional are directly adjusted to 
reproduce at best the nuclear properties. The direct fit to experiment demonstrates 
that the functional contains already at the so-called mean-field level complex many- 
body correlations. The merging of terminology like HF (TDHF) and EDF (TD-EDF) 
stems from a useful artifact of the nuclear mean-field where, most often, the EDE 
is constructed using the concept of effective interaction. It should however be kept 
in mind that the effective hamiltonian technique is actually just a practical way to 
construct an energy that becomes a functional of the density. The difference between 
EDE approach and the Hamiltonian case becomes evident noting that the effective 
interaction vertices might themselves depend on the density of the state to which the 
effective hamiltonian is applied or noting that the interaction in the mean-field chan¬ 
nel differs from the one in the pairing channel. This subtlety has always to be kept 
in mind when considering mean-field or beyond mean-field approaches for nuclei. Eor 
instance, it is safer to directly formulate the variational principle starting directly from 
the energy E{p) instead of the hamiltonian H. If pairing is neglected, one possible 
action to minimize would be: 

r^f 

S= dt 

Jti 

where the single-particle state \ipot) are the occupied state, i.e. p{t) = 

Then, the mean-field is directly defined in terms of the functional derivative of the 
energy h[p\ = dE[p]/dp. 




( 10 ) 


2.1 Macroscopic reduction of the mean-field dynamics 

The starting point of macroscopic approaches to small and large amplitude collective 
motion is to select few collective variables that are anticipated to play a dominant role 
in the physical phenomena under interest. Eor small amplitude vibrations, the operators 
generally correspond to multipole deformation operator [29]. When more than one nucleus 
is concerned like in fusion and/or transfer reactions or larger amplitude deformation, a 
variety of collective variables are introduced: relative distance between nuclei, elongation, 
neck, mass/charge asymmetry... Eor simplicity, a single collective variable q is considered 
here. In a simple markovian limit, the classical evolution of a macroscopic DOE might be 
schematically given by: 


dq p 

dp fdV{q)\ lfdM\, 


( 11 ) 
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Figure 1: Schematic illustration of the adiabatic (left) and diabatic (right) motions. In the 
adiabatic limit, in case of single-particle crossing, the lowest levels are always occupied. In 
the diabatic limit, the same levels remains occupied independently of level-crossing. This 
is for instance what happen if a simple scaling assumption is made for the wave-function. 
This scaling assumption is usually the starting point of hydro dynamical description. 


In these equations, p is the conjugate variable to q. M{q) and V{q) are respectively the 
collective inertia and potential energy landscape. The dissipation term j{q).q is generally 
introduced to account approximately for the coupling between q and other internal degrees 
of freedom, while 5p is a fluctuating term consistently introduced to fulfill the fluctuation- 
dissipation theorem requirement. In particular it might include possible thermal fluctua¬ 
tions. While simple approximations can be made, large effort is made to put as much as 
possible microscopic aspects in the different ingredients entering into the macroscopic equa¬ 
tions of motion |30l|3l] . There are two main difficulties in the macroscopic framework (i) 
one should a priori pre-select very few collective DOFs and therefore, one might miss im¬ 
portant effects that are not anticipated. The validation of the hypothesis is generally made 
a posteriori by confronting results with experiments, (ii) The link between microscopic and 
macroscopic description can be made in general assuming specific approximation for the 
microscopic case, like very slow collective motion (adiabatic limit), or completely diabatic. 
An illustration of the two limits in the case of level crossing is given in Fig. 

One of the great advantage of dynamical mean-field approach is to avoid arbitrariness 
in both (i) and (ii). Indeed, it does not pre-suppose that any one-body DOF plays a domi¬ 
nant role with respect to the others and do not assume a priori adiabaticity or diabaticity. 
Interesting connections between TD-EDF and the classical collective dynamics given by Eq. 
0 can be made by considering specific limit of the evolution, either slow collective motion 
and/or hydro dynamical evolution. The starting point of the adiabatic limit to TDHE is to 
assume that the one-body density takes the form [32j |33] 

p{t) = (12) 

where both x(t) and po(^) are hermitian matrices. x(t) plays the role of a velocity field that 
is common to all single-particle states. In the case of slow motion, x(t) is small and the 
exponential can be expanded. The variational principle § then transforms as: 

S=[ dt{-l^[po{t)x{t)]-E[po{t),x{t)]} , (13) 

Jto 

that is similar to the action in classical mechanics. The connection with the macroscopic 
picture is then achieved by assuming that the density p(t) identifies with the density along 
a potential energy landscape obtained by minimizing the EDE constraining few collective 
DOE go,, i.e. E{po{t)) = E{{qa{t)}). x is then directly proportional to the velocities q^. 
The Adiabatic limit of TD-EDE provides practical formula to evaluate the collective inertia 
and/or the collective potential. It however misses completely dissipative aspects since only 
few DOEs are explicitly followed in time. Possible coupling, even with other one-body DOEs, 
are therefore automatically neglected. 




Beyond Mean-Field 



Figure 2: A schematic representation of the different space of observables considered in 
the present article. The red arrows indicate the coupling between different spaces that will 
induce dissipation. Within the space of one-body DOF, focusing on few collective DOF, we 
expect a coupling between the collective variables and other one-body DOFs. This will lead 
to the one-body dissipation mechanism that is already contained in microscopic mean-field 
theory. Going beyond mean-field can be considered as including two-body or higher-order 
correlations and/or quantum fluctuations in collective space. Treating the coupling between 
one-body DOFs and new DOFs is at the heart of dissipation phenomena induced by internal 
correlations. 


The full TD-EDF includes this coupling and therefore should contains a priori the one- 
body dissipation mechanism (see the schematic picture given in Fig. |^. However, in the 
general case, the microscopic dynamic seems difficult to connect analytically to the transport 
coefficients of Eqs. (11). Nevertheless, a method was developed in Refs. [TTJ [S!]. to 


extract directly the mass parameter, collective potential and dissipation coefficient for the 
relative distance during fusion. The method can be described as follow. Given any one- 
body observable Qa- Its dynamical evolution can directly be obtained along the dynamical 
path from qa{t) = TT{Qap{t)). Eor the relative motion, both the relative distance R{t) and 
associated relative momentum P(t) are known at all times. Assuming that their motion 
match the classical Eqs. 0- one can invert these equations to get the mass, the potential 
and the dissipation kernel. Note that, since TD-EDE is a deterministic approach, there is 
no fluctuating contribution and Spa = 0. The very reasonable collective potential extracted 
with this method and the compatibility of dissipation rate compared to empirical approaches 
again point out that the mean-field dynamic leads to a quasi-classical evolution in collective 
space. 

TDHE and/or TD-EDE is not expected to carry the information on two-body observ¬ 
ables. Indeed, the theory is built up as a functional of the one-body density and its aim 
is solely to provide an accurate description of this quantity. Mean-field is therefore not 
built-up to carry the information on the fluctuations in collective space since it belongs to 
two-body observables. Indeed, taking two one-body observables Qi and Q 2 , their quantum 
correlations are given by: 


^QiQ2(^) ~ -^{QiQ2 Q2Q1) - {Qi){Q2), 


(14) 


that obviously involves the two-body operators Q1Q2 and Q2Qi^ In a mean-field framework, 
these fluctuations become also a functional of the one-body density matrix and can be 
expressed as: 


= Tr[QiK*)Q2(l-p(i))] 


( 15 ) 
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Figure 3: Illustration of the energy sharing during the fusion of ^^0+ assuming a 
simple macroscopic reduction of the mean-field dynamic. Adapted from Refs. [m |34] , the 
dissipated energy reported here is given by the simple formula F^diss = Jq 


where the factor (1 — p) reflects the Fermi statistics of nucleons. It is worth mentioning that 
this fluctuations are a priori non-zero even at zero temperature. Therefore, the mean-field 
theory cannot completely be considered as a fully classical theory. The situation we are facing 
with mean-field is similar to the case of coherent states in quantum mechanics. Coherent 
states are gaussian wave-packets whose properties is fully determined by the knowledge 
of their mean position q and mean momentum p. Using coherent state in a variational 
principle leads to classical equation of motion for {q{t)^p{t)). However, even if the quantum 
fluctuations are not completely neglected, most applications of mean-field theory to nuclear 
dynamics points out that these fluctuations are severely underestimated compared to reality. 
The treatment of correlations and/or fluctuations requires to go beyond the independent 
particle or quasi-particle picture [35]. 

3 Beyond mean-field approach to large amplitude col¬ 
lective motion 

Beyond mean-field approaches is usually facing the difficulty to be much more involved 
in practice than the mean-field level. A common idea to BMF theory is to generalize mean- 
field by including more complex internal degrees of freedom linked to two-body or higher 
correlations. If we denote generically by Q = {(Qa)} the expectation values of the set of 
one-body DOFs and by B = {{Bi)} other degrees of freedom, the exact equation of motion 
can be formally written as: 

[ ^(Qa) = {[Qa,H]) = 


= {[Bi,H]) = GiiQityMt)) 


(16) 
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Figure 4: Schematic illustration of the different scenarios to include effect BMF keeping only 
the one-body density as relevant DOF. Starting from the deterministic mean-field evolution, 
different degrees of sophistication can be used to include beyond mean-field effects. Keeping 
only the mean-field and the dissipative term V in Eq. (19) lead to a new deterministic 
picture where the effect of complex DOFs is treated in an average way. Keeping both the 
dissipation and fluctuating term 5J^ leads to a quantum Boltzmann approach where both 
dissipation and fluctuations are introduced consistently (bottom). A similar picture was 
used in ref. [36] to describe semi-classical methods dedicated to heavy-ion collisions around 
Fermi energy. 









The main difficulty in solving this problem is the catastrophic explosion of number of DOF 
to follow in time. Indeed, supposing that the number of one-body DOF is the number 
of correlation matrix elements is N‘^I2, the number of third moments is /2^ ... The set 
of equations given above contains the coupling between one-body DOF and more complex 
DOF that is responsible of two-body dissipation (see schematic picture |^. 

Except very rare cases, specific strategies should be used to reduce the complexity of this 
problem. The great simplification of mean-field approximation is that all DOFs B becomes 
functional of the relevant DOFs Q, i.e. (B) = /((Q)). Then, the second set of equations in 


(16) becomes redundant and we simply have N coupled transport equations: 


This equation is nothing but Eq. already shown in previous section. 

Three strategies commonly used to extend mean-field are summarized below: 


(17) 


(i) Eollowing all DOEs in the coupled equations (16) is usually prohibitive. When possible. 


it is preferable to not follow them explicitly but correct the {Qa) evolution by including 
the effect of other DOEs in an approximate way. The problem we face is similar to 
an open quantum system (OQS) problem where a system is coupled to a complex 
environment m This similarity provides strong guidance to the N-body problem. 
The starting point to avoid the explicit treatment of the environment DOEs is to first 


integrate the second equation in (16), leading to: 


= (Mto)) + / &(Q(s);B(s))ds 


/■ 

Jto 


(18) 


The appearance of the expectation value of B at different time on the left and right 
hand side of this equation is similar to what happen in the Green function theory. 
If the coupling between different degrees of freedom is not too strong, a perturbative 


approach can be constructed from Eq. (18). In many-body physics, an expansion a 
la Nakajima-Zwanzig (NZ) [39l [40] is followed. Reporting the integral equation (18) 
with some approximation in the first equation (16), lead to new equation of motion 


for the relevant DOE only, that could be written as: 

^=T{Q) + V{Q) + 5T{Q). 


(19) 


The first term in this equation is nothing but the mean-field. The second term, some¬ 
times called dissipation kernel or dissipator contains the accumulated effect of the 
coupling from the initial time to to a given time t. It stems from the integral in time 
in Eq. (18) and therefore depends on the full history of the system, i.e. on Q(s) 
with to < s < t. This non-locality in time is sometimes called memory effect and/or 
non-Markovian effect. The importance of memory is well known especially in quantum 
systems [37] and has been for instance illustrated in Ref. [41]. It is worth mentioning 
that the Nakajima-Zwanzig technique that leads to equations that are non-local in 
time for the relevant DOE has been recognized to not be the most efficient way to 
treat open quantum systems. Eor instance, other techniques like time-convolutionless 
(TCL) seems more adequate [37]. In that cases, the equation of motion for relevant 
are time local and memory effects are incorporated in the transport coefficients. This 
issue is largely unexplored in our field. 


The analogy between the N-body problem and OQS can be further exploited to obtain 
a Lindblad equation from the dissipation kernel starting from Eq. ( [l^ [38] . 
Ultimately, the dissipation kernel can then be interpreted in terms of quantum jumps 
between independent particle states. 


The third term contains the effect of initial correlations that are propagated in 

time through the mean-field. Invoking the absence of knowledge of this initial corre¬ 
lation this term is sometimes treated as a noise added to the mean-field. Comparing 










the original mean-field equation 0 and the new one, Eq. ([T^ that includes even¬ 
tually dissipation and fluctuation, the increase of difficulty stems from the memory 
effect that is important in quantum mechanics and from the necessity to follow several 
trajectories instead of one if is treated as a fluctuating term. The different 

scenarios of beyond mean-field dynamics without treating complex degrees of freedom 
explicitly are shown in Fig. 

(ii) Another possibility is to add new degrees of freedom and try to solve a generalization 
of mean-field equation in the form (16). Due to the number of DOFs (one-body, 
two-body, three-body, ...) it is just impossible to include all of them and a specific 
truncation is required, for instance in the Time-Dependent Density Matrix approach 
(TDDM) [5], the two-body correlation matrix is followed in time while neglecting 
three-body correlations. The absence of a clear criteria to truncate the hierarchy of 
equations jeopardize applications (see for instance Ref. m)- Still few applications 
have shown that correlations can be partially treated in this way (see 03111 ] and 
references therein). 


(iii) Finally, a third possibility is to consider more complex wave-functions allowing for a 
mixing of single-particle states. For instance, if the aim is to re-quantize the motion in 
collective space. A possible method is to use trial state vector that is already a mixing 
of several independent particle or quasi-particle states. Then, the trial wave-function 
to be used in the variational principle can generically be written as: 




( 20 ) 


where {|^(Q, t))} stands for a set of eventually time-dependent quasi-particle vacuum. 
This general approach is called Time-Dependent Generator Coordinate Method (TD- 
GCM). While this approach starts to be a standard tool for nuclear structure, very 
few discussions and applications have been made for the dynamical process (see for 
instance the discussion in ref. [l^ |46] and application in [47]). In particular, one of 
the extra difficulty is the doubling of collective space size compared to the static case, 
due to the necessity to include pairs of conjugated collective variables. 

The two strategies (i) and (ii) have attracted more attention in the last 30 years in the 
non-equilibrium context compared to the strategy (iii). However, strategy (iii) is the one 
that is now being recognize as one of the standard tool to extend mean-field and obtain 
spectroscopic information in nuclei. In the latter case, important discussions have been 
made recently to reconcile the GCM approach that is well defined starting from a true 
Hamiltonian [3| with the energy density functional theory where no Hamiltonian exists at 
the first place iMiiiaiioiini]. 

The TD-EDF theory is also strongly guided by the Hamiltonian case. Most often, if 
not always, final equations of motion in the three strategies are derived starting from a 
Hamiltonian and directly applied within the EDF context. The analogy with the Hamil¬ 
tonian case should ultimately be taken with care and in particular, strategies (i), (iii) and 
(iii), when they are used in the nuclear context, should be seen as methods to enrich the 
functionals. Indeed, in the TD-DFT context, the very notion of mean-field level and beyond 
mean-field level is a ill defined concept since already the approach should provide an exact 
description of the local density without invoking additional DOFs. In the DFT context, the 
three strategies listed above should be regarded as follow: 

• The strategy (i) can be seen as a way to construct a richer functional of the one-body 
density matrix. In particular, the appearance of time non-locality is not surprising in 
the DFT context where the functional is expected to depend on the history of the local 
density [52]. If the non-locality is important, the functional is called ”non-adiabatic”. 
This terminology should not be confused with the non-adiabaticity discussed in section 

\Tn 




• The strategy (ii) should be interpreted as an alternative functional theory where the 
energy becomes a functional not only of the one-body density but also of more complex 
degrees of freedom. 

• The strategy (hi) is certainly the one that is more difficult to understand in terms of 
a DFT. In Refs. [53l [54], it has been shown that the many-body state |T) entering 
in Eq. ( [^ can be regarded as a specific trial state from which the one-body and 
two-body densities are constructed leading to an energy functional of these quantities. 
We then end up with the schematic DFT scheme: 

1^) —^ (p, P12,-")— > pi2,' ")• (21) 

The work in [ssiisi] was dedicated to static properties and the application to nuclear 
dynamics within the EDF approach still remains to be formally clarified. 

In this section, a schematic presentation of beyond mean-field approaches is given trying 
to underlined some specific discussions related to the use of the Energy Density Functional 
approach. As we have seen above, the method used to extend mean-field is not unique and 
depends essentially on the physical process we want to include. Below, I concentrate on the 
specific description of finite many-body systems at low internal excitation. In that case, it is 
expected that the main difficulty is to treat quantum fluctuations beyond the independent 
particle approximation. Discussions on methods specifically dedicated to system at higher 
internal excitation can be found in [5l |211 [55] . 

4 Propagation of Quantal fluctuations in Beyond Mean- 
Field theories 

One of the deficiencies of mean-field transport theory is its inability to describe quantal 
fluctuations around the average properties in collective space. This shortcoming is directly 
related to the quasi-classical nature of the collective evolution. At low internal excitation, 
this aspect is anticipated to be the most important effect beyond mean-field. The treat¬ 
ment of quantum fluctuations requires a priori to re-quantize the collective dynamics by 
considering a quantum wave-function in collective space, like in the TD-GCM (item (iii) 
above). This approach, involving a coherent mixing of many independent quasi-particles 
many-body states remains quite involved and has rarely been exploited in the dynamical 
context. Alternative methods have been proposed that are trying to account for quantum 
fluctuations while keeping the simplicity of independent particle/quasi-particle states. Some 
of the approaches, depicted schematically in Fig. [^are briefly described in this section. 

4.1 The Balian and Veneroni variational principle: 

The variational principle proposed in the 80^^ by Balian and Veneroni, Eq. Q can be 
used to derive time-dependent mean-field if the relevant degrees of freedom are restricted to 
one-body DOFs. Generalization of TDHF can be obtained either by considering additional 
DOF, like two-body or higher correlations or by considering more general many-body density 
matrix or both. There are essentially two difficulties in using the BV variational principle (i) 
if complex density matrix are used, one looses the advantages of the the algebra of densities 
that are given by Eq. <11 , see for instance m (ii) the equations of motion obtained by 
adding more complex internal DOFs become rapidly rather involved [24l[25j[26] and, as far 
as I know, have not been so far really integrated in time. In order to include the fluctuations 
associated to a given one-body observable without facing the problem (i), BV proposed to 
optimize the set of observables given by A = while D is still given by Eq. Then, the 
expectation of A becomes the generating function of all moments of Q. Still, the equations to 
be solved remain complicated. However, it was shown in Ref. [25][56] that, in the specific case 
where the density identifies with a Slater determinant density, the fluctuations associated to 
a given variable can be obtained by performing uniquely forward and backward mean-field 
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Figure 5: Schematic illustration of the different approaches discussed in the text to describe 
quantum fluctuations beyond the mean-field theory (see text for more details). 






evolution between initial and final times ti to tf. The forward evolution is obtained from an 
initial one-body density p{ti) and leads to the density p{tf). Then a second evolution starts 
from a density r]{tf) = p{tf)e~^^^ that is evolved backward to give p{ti) (see illustration 

in Fig. 1^. Finally, the fluctuations at final time tf that accounts partially for two-body 
effects are given by: 


4 «/) = ( 22 ) 

This approach can eventually be generalized to include more than one collective variable 
and also applies if the state has non-zero pairing. With the progress made to perform TD- 
EDF, a renewal of interest is observed to extend the mean-field through the BV variational 
principle [58l|59l|60] (see also the contribution of C. Simenel in this book). In particular, it 
is expected that, in a regime where the backward propagation remains close to the forward 
evolution, the quantum fluctuations are corrected due to the Random-Phase-Approximation 
(RPA) correlation around the mean-field trajectory. However, there are situations where a 
small difference in the initial conditions is anticipated to lead to large deviations even over 
short time. This is the case, when the system is close to a bifurcation point associated to 
a spontaneous symmetry breaking. While this issue remains to be clarified, the application 
shown in ref. m seems to indicate that the BV variational principle is inadequate in that 
case. 

4.2 The time-dependent RPA theory 

A closely related approach that has been proposed approximately at the same time as 
the BV method consists in building up directly a RPA like theory all along the mean-field 
trajectory [62j|63]. This is the so-called time-dependent RPA theory. In the RPA (or QRPA) 
approach, a correlated ground state is constructed starting from a static independent/quasi¬ 
particle vacuum. This correlated state becomes then a vacuum for quasi-boson operators 
written as a coherent superposition of particle-hole excitations [3]. In RPA, the collective 
energy landscape is assumed to be a harmonic and quantal fluctuations beyond mean-field 
can be extracted. In Ref. [62], this approach has been generalized to include quantal 
fluctuations along the mean-field trajectory. Assuming that (i) the average trajectory is 
weakly affected by correlations and identifies with the TDHF trajectory (ii) that the quantal 
fluctuations remains small, a RPA theory is constructed at each time. In this theory, a new 
state |T(t)) is constructed that is a vacuum of a set of quasi-boson operators {Ba{t)}, i.e. 
Ba{t)\'^{t)) = 0. To make contact with the collective space evolution, it is convenient to 
write the quasi-boson operators as: 

Bait) = Bait) - iXaQait) (23) 

where = “^iPa)- The operators iPa,Qa) are conjugated operators along the path. In 
particular, they defined specific trajectories around the average mean-field evolution through 
the transformation: 


|[$ + 5$](t)) = (24) 

where are small parameters. 

The equation of motion are obtained using the variational principle for the correlated 
state |T(t)) leading to coupled equations between the state |^(t)) that evolves through 
the standard mean-field assuming (e^ = 0, = 0) and the set of operators evolution 

{Pa{t),Qa{t)), written in the Heisenberg representation. Technical details are given in ref. 
[62] . As in the RPA case, it could be shown that the state |T(t)) is a coherent mixing 
of two-particles two-holes excitation with respect to the independent particle state |^(t)) 
and therefore contains additional correlations leading to enhanced collective fluctuations. 
Again, as in RPA, the correlated state do not need to be explicitly constructed thanks to 
the quasi-boson approximation [3]. 



It is finally interesting to note that some connection can be made with other approaches. 
Indeed, let us assume a practical situation where a set of DOF {Qa{ti)^ Pa{tf)) are known 
at initial time, as well as the TDHF starting point \^{tf)). At all time, the operators 
{Qa{t), Pa{t)) are defined through the implicit formulas: 


Pa{t)\m) 

Qc.{t)\m) 


i lim 
—i lim 


6a 




(25) 


where \^Sa, (^)) and (t)) are independent particle/quasiparticles state that evolves through 
their own mean-field with an initial condition that is slightly shifted either in or Qq,. This 
method is illustrated in Fig. Again this method only applies if the deviation between the 
different states considered initially remains small during the time evolution. 


4.3 The stochastic mean-field theory 

One of the advantage of the two approaches depicted above is that they only requires 
to follow standard mean-field. The practical way that consists in following two mean-field 
trajectories that are close in collective space is anticipated to fail when large deviation will 
occur even if the initial conditions are very close from each other. This is the case for instance 
in quantum chaotic systems or close to a bifurcation point associated to the appearance of 
saddle points in collective space. Another important aspect, is the inability of mean-field 
to spontaneously break symmetries. Such symmetry breaking is of particular importance 
especially in quantum many-body systems and induces usually large quantum fluctuations 
in collective space. 

In recent years, a new approach [Ml [21], called Stochastic Mean-Field (SMF) has been 
proposed to overcome these difficulties where also only mean-field evolution is needed. The 
two main hypothesis of SMF are: 

1. The mean-field equation of motion of the density can be regarded as the classical or 
quasi-classical limit of a many-body fermionic problem. As we have seen above, mean- 
field equation leads indeed to quasi-classical trajectories where many quantum effects 
are missing. 

2. Quantum dynamic and in particular quantum fluctuations evolution can be simulated 
using classical trajectories with a proper sampling of the initial collective phase-space. 
Note that such a classical mapping is a known technique to simulate quantum objects 
and might even be exact in some cases [65] [66]. While a different method is used 
to treat the phase-space, a similar strategy has been used in bosonic systems, the 
so-called truncated Wigner method m- 

In practice, SMF has been directly formulated in density matrix space. Let us assume 
that the aim is to improve the description of a system that, at the mean-field level and time 
to, is described by a density of the form: 

Pi^o) = 'p \‘Pi{to))ni{(pi{to)\. (26) 

i 

Then the mean-field evolution reduces to the evolution of the set of single-particle states 

= h[p]\<fi{t)), (27) 

while keeping the occupation numbers constant. 

In the SMF approach, a set of initial one-body densities 

P^po) = 


( 28 ) 




is considered, where A denotes a given initial density. The density matrix components 
are chosen in such a way that initially, the density obtained by averaging over different initial 
conditions identifies the density (26). 

It was shown in ref. [64] that a convenient choice for the statistical properties of the 
initial sampling is 


ptj{to) = SijJii + Sp^j{to), 

where 6pfj{to) are mean-zero random Gaussian numbers while 


(29) 


(30) 


The average is taken here on initial conditions. In this approach, each initial condition given 


by Eq. (28) is then separately evolved with its own mean-field: 


ih 


dp^{t) 

dt 




(31) 


independently from the other trajectories. Therefore, the evolution along each trajectory 
is similar to standard mean-field propagation and can be implement with existing codes. 
Two schematic illustrations of the stochastic mean-field method are given in bottom parts 
of figure]^ Note that the SMF approach has been generalize to include pairing in ref. [68] . 
The technical details and success of SMF approach have been recently reviewed in ref. m 
and only specific aspects related to collective motion are given here. 


4.3.1 Collective dynamic in SMF 

In the original mean-field approach, even if the theory is underestimating the quantum 
fiuctuations compared to reality, the mean values and fiuctuations of a given one-body 
collective observable O are still given by a quantum average, i.e. 

(Q) = Tr{Qp), al = {Q^)-{Qf =Tr{QpQil-p)). (32) 

In SMF, along a given trajectory associated to p^{t)^ the observable evolution is given by: 

Q\t) = T,{Qp^{t)). (33) 


The very notion of quantum average is lost and the variable Q^{t) is now considered as 
a classical object. For instance, fiuctuations are estimated using the classical average over 
different trajectories, i.e. : 

= W-Q‘' (34) 


where A^traj is the number of trajectories. It is noteworthy that the sampling of initial 
conditions made under the two constraints given by Eqs. (29) and (30) insures that the 
classical average Eg made over the initial conditions on SMF match the quantum average 
cTg of the original quantum mean-field at initial time. 


4.3.2 SMF with small deviation around the mean-field evolution 

One direct proof that a classical mapping of mean-field theory with appropriate sam¬ 
pling of initial conditions can grasp effects beyond the independent particle or quasi-particle 
picture, is the possibility to make connection with the BV and/or TDRPA when trajectories 
remain close to the average mean-field evolution all along the trajectory (third panel case 
in Fig. [^. In Ref. [26|, it has been shown that, assuming small deviation with respect to 










the mean-field evolution, the BV variational principle corrects the quantum fluctuations by 
including the TD-RPA correction. Then the quantum fluctuations at final time are given 
by: 




= Tr [Q{ti)p{ti)Q{ti){l - p{ti))\ 


(35) 


where Q{ti) is the observable that is propagated backward from tf to U through the time- 
dependent RPA equation with initial condition Q{tf) = Q. In [64], assuming small fluc¬ 
tuations, it was also shown that the SMF theory leads also to Eq. (35) showing that it 


can describe at least TDRPA correlations beyond mean-field. This was verified explicitly 
for the Lipkin model case [69] when the collective energy landscape is harmonic. In that 
case, both the BV/TDRPA theory [61] and SMF gives identical results that can hardly be 
distinguished from the exact solution. 

The implementation of SMF is however a priori more demanding than the practical 
solution of the BV approach described in section |4.1| or |4.2| since it requires to consider 
enough trajectories to reproduce the initial phase-space. However, similarly to the BV or 
TDRPA, assuming that all trajectories remains close to the average mean-field (TDHF) 
evolution, the explicit sampling of trajectories can be avoided. Then, the numerical cost to 
obtain corrections to mean-field within SMF becomes of the same order as for the original 
mean-field. The simplified implementation of SMF can be summarized as follow. Starting 
from the mean-field equation, written in a generic way as in Eq. 0: 


dt 




(36) 


and focusing on a specific variable Q, following the same strategy as for TDHF (Eq. 0) 
a macroscopic reduction of SMF can be made leading to an equation of motion: 


dQ^ 

dt 




(37) 


where F is the driving force associated to the collective potential, z/ is the dissipative kernel 
that contains on-body dissipation effect. Both F and u are already present in the mean-field 
level. 6Q^{t) is the additional fluctuating term that stems from the initial sampling. It is a 
priori non-Markovian. Assuming that all trajectories remain close to the average trajectory 
and that this trajectory is identical to the mean-field one, one can eventually approximate: 


F{Q\t)) ~ F{{Q)), (38) 


where (Q) is the expectation value associated to the standard time-dependent evolution. 
With this simplification, the average evolution of the fluctuation in SMF ^Q{t) contains an 
additional diffusion term associated to the classical fluctuations For a given 

observable, the estimate of this quantity is far from being straightforward. However, using 
the Wigner transform and assuming no memory effect (markovian limit), an approximate 
method has been proposed in Ref. HO] and further applied in [zn na ESI to heavy-ion 
collisions. The great interest of this approach is that corrections beyond mean-field are 
obtained directly from the mean-field evolution. Very recently, it was shown that a fully 
quantal treatment including memory effect can be made m opening new perspectives. 


4.3.3 SMF with large deviation around the average 

The above simplified method will not work if the densities along each trajectory starts 
to strongly deviate from the average. An illustration of such a case is given in the bottom 
part of figure [^ Then, a direct sampling of the initial phase-space followed by a set of time- 
dependent mean-field trajectory is unavoidable. Up to now, the SMF theory in its full glory 
has been applied to rather simple models. The application to the Lipkin model [69] was the 
first attempts and has shown that not only the SMF theory is equivalent to the TDRPA for 
harmonic motion around the ground state but also that it works in situations where other 






approaches fails, for instance to describe the non-equilibrium dynamic close to a spontaneous 
symmetry breaking associated to a quantum phase-transition. In ref. [68], this aspects has 
been used specifically to take advantage of the U{1) symmetry breaking associated to particle 
number conservation, leading to a generalization of SMF including pairing. The new theory 
was validated using the pairing hamiltonian. The capability of SMF has then been further 
demonstrated with the Hubbard model, where fast non-equilibrium process can occur after 
a quantum quench. In that case, the approach was competitive with state of the art non- 
equilibrium Green Function theory m- The successes of SMF are very promising. Still 
many aspects needs to be clarified: 

• Predictive power: In different applications of SMF, it was systematically observed 
that this theory can provide a very good approximation beyond mean-field in the weak 
coupling regime, i. e. when the residual interaction is weak. In some cases [6^ [68] . 
the SMF results cannot be distinguished from the exact solution. When the strength 
of the two-body interaction increases, usually the SMF approach is rather good over 
a certain time and then, starts to deviate from the exact evolution. In the Hubbard 
model case, in the strong interaction regime, it was observed that the SMF fails to 
reproduce the dynamics even at rather short time scale. Up to now, the approach has 
been tested in a rather empirical way and its success can only be tested comparing 
with the exact solution. In the near future, it would be desirable to seek for some 
criteria that would allow to anticipate when the theory would have a good predictive 
power. 

• Symmetry breaking: A clear interesting aspect of the approach is to allow for initial 
densities that can break some symmetries that are respected by the initial average 
density p{t). In a physical system many symmetries can be broken. In practice, the 
density p^ can be chosen to break all possible symmetries including U{1) symmetry or 
can break only selected symmetries while some others are respected. As a consequence, 
the SMF theory is not unique and can be implement at various levels of symmetry 
breaking. Obviously, the more symmetry are broken, the more demanding is the 
calculation. The results of SMF strongly depend on the symmetries that are initially 
broken. For instance, in the Hubbard model, a better result is obtain if the spin 
up/spin down symmetry is broken compared to the case where it is not broken. On 
opposite, no improvement was observed if the U{1) was broken compared to the case 
where it was not. Again, the choice of breaking or not of some symmetries is rather 
empirical and requires a priori some physical intuition. 

• Gaussian sampling of the phase-space: In the original SMF theory formulation, 
only the first and second moments of the initial density fluctuations are constrained to 
match with the initial quantum fluctuations. This hypothesis, that implicitly assume 
a Gaussian initial distribution of the density matrix has been recently questioned in an 
exploratory study Hg. Using a more realistic sampling of the initial phase-space, it was 
shown that the SMF dynamics can be further improved, especially for larger coupling 
and longer time. This clearly points out that the approach and its predictive power is 
largely unexplored. A key challenge to further progress is to have an efficient method 
able to provide a realistic phase-space in complex many-body fermionic systems. 

• Correlated initial state: Again in [76|, it was also shown for the first time that the 
SMF can be applied even if the initial state is not an independent particle system but 
contains correlations. The key of the success was the realistic sampling of the initial 
phase-space. 


5 Summary 

After the first application in nuclear physics of time-dependent mean-field based on 
Skyrme effective interaction, it was realized very rapidly that this approach should be ex¬ 
tended to include beyond mean-field effects. During many years after this first application, a 


great number of famous scientists have proposed a variety of approaches that includes differ¬ 
ent effects [35] . Most of the frameworks turn out to be so complex that they have never been 
applied so far. One of the difficulty is the quantum nature of nuclei. In view of this com¬ 
plexity, some theories where applied using a semi-classical approximation [saizzi- In recent 
years, great progress have been made in the application of TD-EDF to the nuclear many 
body problem. This has renewed also the interest of developing quantum transport theories 
beyond mean-field. In the present notes, an overview of several techniques that incorporate 
some effects beyond the independent particle/quasi-particle approaches is given. 

In the second part of this article, a focus is made in finite many-body systems at low 
internal excitations. In that case, quantum fluctuations in collective space are expected 
to play an important role. Mean-field approaches based on one-body density matrix are 
known to poorly describe quantum collective fluctuations. Theories that can be regarded as 
practical tools to improve mean-field in the context of large amplitude collective motion are 
described: namely the Balian-Veneroni variational principle, the Time-Dependent Random 
Phase Approximation and the Stochastic Mean-Field theory. Their advantages and draw¬ 
backs are discussed. Among them, the Stochastic Mean-Field approach seems to provide a 
versatile tool able to describe both the collective motion in the harmonic limit as well as 
strong non-equilibrium process occurring close to a bifurcation point or a quantum phase- 
transition. Last, I would like to mention that the present discussion is by no way exhaustive 
and the future of transport theories applied to finite quantum many-body systems is largely 
open to new ideas. 
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